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In this work we study how scalar fields may affect solar observables, and use the 
constraint on the Sun's central temperature to extract bounds on the parameters 
of relevant models. Also, a scalar field driven by a suitable potential is shown to 
produce an anomalous acceleration similar to the one found in the Pioneer anomaly. 



1. Introduction 

Scalar fields play a fundamental role in contemporary physics, with ap- 
plications ranging from particle physics to cosmology and condensed mat- 
ter. In cosmology, these include inflation 1 , vacuum energy evolving and 
quintessence models, the Chaplygin gas dark energy-dark matter unifica- 
tion model and dark matter candidates (for a complete list, see [2] and 
references within). A scalar field has recently been suggested as a possi- 
ble explanation for the anomalous acceleration measured by the Pioneer 
spacecraft 3 . On an astrophysical context, a scalar field could be present as, 
for example, a mediating boson of an hypothetical fifth force. This can be 
modelled by a Yukawa potential, given by Vy(r) — Ae~ mr /r, where A is 
the coupling strength and m is the mass of the field, setting the range of 
the interaction Ay = mT 1 . The bounds on parameters A and Ay = mT 1 
can be found in Ref. and references therein. In this work, a scalar field 
with an appropriate potential is shown to account for the anomalous Pio- 
neer acceleration. Also, the effects of this and other scalar field models on 
stellar equilibrium is reported. 
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2. Scalar field induced Pioneer anomaly 

The hypothesis that the Pioneer anomaly is due to new physics has been 
steadily gaining momentum. Fundamentally, this arises because is seems 
unfeasible that an engineering explanation, affecting spacecrafts with dif- 
ferent trajectories and designs, could lead to similar effects. The lack of 
a simple explanation could hint that this anomalous acceleration is the 
manifestation of a new force. 

Given the broad range of applicability of theories with scalar fields, it 
is possible that such an entity is responsible for the anomaly 3 . A natural 
candidate would be the radion, a scalar field present due to oscillations of 
inter-branc distance in braneworld models. However, this does not account 
for the phenomena 3 . Fortunately, one can put down a model with a scalar 
field that doesn't stray away much from the already accepted theories aris- 
ing from cosmology. Specifically, we assume a scalar field with dynamics 
ruled by a potential of the form V{<j>) = —A 2 cf)~ 2 , where A is a constant. 
The scalar field obeys the equation of motion 

4>"{r) + V(r) = 2A 2 ^ 3 , (1) 
and admits the solution 



0(r)= (jj s/Tr = (3-^ . (2) 
The — linearized Einstein's equations reads 

/"(r) + ^/'(r) = 2K(l-^)^ 2 /3 2 r- 1 (3) 



whose solution is 



/M'Vl^-^Wry • < 4 » 

The resulting acceleration felt by a test body is given by 



C jSAn jSACn 
T 2 ~ V 2T + V 2~ 
The first term is the Newtonian contribution, and identifying the 
second with the anomalous acceleration cla — 8.5 x 10~ 10 ms~ 2 , sets 



a r = -~2 - \l o— +Vo— • ( 5 ) 
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A = 3.26 O/i/k = 4.7 x 10 42 m~ 3 . The last term is much smaller than the 
anomalous acceleration for 4C/r <C 1, that is, for r ^> 6 km; it is also much 
smaller than the Newtonian acceleration for r <C 2.9 x 10 22 km w 100 Mpc. 
It has been shown 3 that the effect of this model on the Doppler effect used 
to calculate the acceleration of the spacecrafts is negligible, indicating that 
the anomaly is real, an not due to a misinterpreted light propagation. 

3. Variable Mass Particle models 

The variable mass particle (VAMP) proposal 5 assumes the presence of yet 
unknown fcrmions coupled to a dark-energy scalar field with dynamics ruled 
by a monotonically decreasing potential. Although this has no minima, 
the coupling of the scalar field <f) to this exotic fermionic matter yields 
an effective potential of the form V e ff{4>) — V{<p) + Xn^cp, where n$ is 
the number density of fermionic VAMPs and A is a Yukawa coupling. In 
the present work we approach a potential of the quintessence-type form, 
V{<p) — u n (f>~ p . As a result, the effective potential acquires a minimum 
and the related vacuum expectation value (vev) is responsible for the mass 
term of the exotic VAMP. Since the number density depends on the 
scale factor a(t), this mass varies on a cosmological timescale. A noticeable 
shortcoming of VAMP models is the presence of weakly or unconstrained 
parameters: the relative density fl^o, the scalar field coupling constant and 
the potential strength. 

One can attempt to overcome this drawback by assuming that all 
fcrmions couple to the quintessence scalar field. This coupling should not 
substitute the usual Higgs coupling, but merely add a small correction to 
the Higgs-mechanism induced mass. Since the effects of the quintessence 
scalar field coupling crucially depend on the particle number density n$, 
one expects this variable mass term to play a more relevant role in a stellar 
environment than in the vacuum. This hints that VAMP model parameters 
could be constrained from stellar physics observables. 

3.1. The polytropic gas stellar model 

The polytropic gas model assumes an equation of state of the form P = 
Kp n+1 / n , where n is the polytropic index, defining intermediate processes 
between the isothermic and adiabatic thermodynamical cases, and K is 
the polytropic constant, which depends on the star's mass M and radius 
R. This model leads to scaling laws for thermodynamical quantities, given 
by p = p c n (O, T = T C 6(Z) and P = P c 9^) n+1 , where p c , T c and P c 
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is the density, temperature and pressure at the center of the star 6 . The 
dimcnsionlcss function depends on the dimcnsionlcss variable £, related 
to the distance to the star's center by r = a£, where a also depends on the 
star's mass M and radius R. The hydrostatic equilibrium condition enables 
a differential equation ruling the behavior of the scaling function #(£), the 
Lane-Emden equation: 



Since the physical characteristics of a star appear only in the definitions 
of the constants K and a, its stability is independent of these quantities, 
and different polytropic indexes n label different types of stars. This scale- 
independence can be related to the homology symmetry enclosed in the 
Lane-Emden equation. The boundary conditions of this differential equa- 
tion are, from the definition p = p c 6(£), 9(0) = 1; also, the hydrostatic 
equilibrium condition implies that IdO/d^]^ = 0. The first solar model 
ever considered corresponds to a polytropic star with n = 3 and was stud- 
ied by Eddington in 1926. Although somewhat incomplete, this simplified 
model gives rise to relevant constraints on the physical quantities. 

The following results are based on the luminosity constraint on the 
Sun's central temperature, AT C /T C < 0.4%. The central temperature can 
be computed from T c = Y n p,GM/ R, where k is the Boltzmann constant 
and Y n depends on the star's mass M and radius R, as well as on the 
dimensionless quantity £i, which is defined through 0(£i) = and signals 
the surface of the star. 

3.2. Results 

Assuming isotropy, a variable mass term leads to both a radial, anomalous 
acceleration = 4>' / 4> < plus a time-dependent drag force ap = — 0/ '<f> < 
2 . The time-dependent component should vary on cosmological timescales, 
and can thus can be absorbed in the usual Higgs mass term. Hence, one 
considers only the perturbation to the Lane-Emden equation given by the 
radial force. We start by defining the dimcnsionlcss scalar field $ = <j>/4>*, 
with (j)* = p crit /2ny, where fiy is the energy density due to the potential 
driving <f> and p cr n ~ 9.48 x 10~ 28 g cmT 3 is the critical density. Assuming 
for simplicity a potential with p = 1, we obtain a perturbed Lane-Emden 
equation 
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where one has defined the dimensionless quantities 



(7) 



GM 



Rc 2 

with W„ = [47r(n 



2.12 x 1CT 6 , C~ l = (n+ l)iV™/( ri+1 )VK„ 1/( " +1) 



d «4 



(8) 



Furthermore, one assumes that the scalar field is only weakly perturbed 
in relation to the cosmological vev of the effective potential. Denoting this 
small "astrophysical" contribution as $<,(£)> one nas ^(0 = ^VM + $ a (£)- 
Hence, the Klein- Gordon equation becomes 



inside the star, and 



2a 2 Xn v Pcq3^ _ 1a 2 \ 2 n v 



Pcrit M 



Pcrit 



(9) 



2a 2 Xn v 



1 - 



2A$ n (fl 



(10) 



Pcrit Pcrit 

in the outer region. In the above, n$ = ny = 3 m~ 3 is the number density 
of fermions in the vacuum. \i is the mean molecular weight of Hydrogen. 
The boundary conditions for the perturbed Lane-Emden are the same as 
in the unperturbed case. For the scalar field, we assume both $(£) and 
its derivative vanish beyond the Solar System (about 10 5 AU). Following 2 , 
one gets that the maximum deviation AT C /T C = 2.82 x 10~ 8 occurs for 
Q,v = 0.7, A = 2.82 x 10~ 14 . Hence, the luminosity constraint is always 
respected as long as one assumes the bound arising from the assumption 
that the "astrophysical" component of the scalar field is much smaller than 
its cosmological vev, A < 10~ 14 . 



4. Yukawa potential induced perturbation 

One now looks at the hydrostatic equilibrium equation with a Yukawa po- 
tential which, after a small algebraic manipulation, implies the perturbed 
Lane-Emden equation 



1 d dO 

edldi 



de 



1 + Ae'^ 1 - jAC n —e-^ 6 
J d£ 



(11) 
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where we have defined the dimensionless quantities 

Cn = V2 7Vy("+ 1 )^ 1 -")/ 2 ("+ 1 ) , 7 = mR (12) 

and £i signals the surface of the star (more accurately, a surface of zero 
temperature, but the difference is negligible. It can be shown that the 
boundary conditions are unaffected by the perturbation 2 . 

One can study the variation of the central temperature as a function 
of A = m" 1 and A, and constraint these so that AT C /T C < 0.4%. The 
parameters were chosen for Yukawa interactions in the range O.li? < Ay < 
10R; the Yukawa coupling A was chosen so that the variation of T c is of the 
same order O(10~ 4 )) as the luminosity constraint. Numerical integration of 
Eq. 11 is then used to derive the exclusion plot of Figure 1, superimposed on 
the accepted bound 4 . Notice that the central temperature is not precisely 
known and it is clear that constraining its uncertainty below 10 would 
yield a larger exclusion region in the parameter space. 

4.1. The Pioneer anomaly 

Following a similar procedure to the one depicted in the above two sections, 
one can prove that a constant, anomalous acceleration a a inside the Sun 
yields a relative deviation of the central temperature which scales linearly 
with cla as 5T C ~ where a© = 274 m.s~ 22 . Thus, the bound 

5T C < 4x 10~ 3 is satisfied for values of this constant anomalous acceleration 
up to auax ~ 10 Og. The reported value is then well within the allowed 
region, and has a negligible impact on the astrophysics of the Sun. 

5. Conclusion 

In this work we have developed a study of the impact of some exotic physics 
models on stellar equilibrium, enabling the extraction of constraints on 
the relevant parameters of the theories. Also, a model exhibiting a scalar 
field with a attractive quintessence-like potential is discussed, which can 
account for the anomalous acceleration felt by the Pioneer 10/11, Ulysses 
and Galileo spacecrafts. 
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Figure 1. Exclusion plot for the relative deviation from unperturbed central temper- 
ature T c , for A ranging from 10~ 3 to 10 _1 , and 7 from 10 _1 to 10 (tip at the top), 
superimposed on the available bounds. 
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